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Abstract 

We calculate critical couplings and string tensions in SU(2) and SU(3) pure lattice gauge 
theory by a simple and inexpensive technique of two-lattice matching of RG block transfor- 
mations. The transformations are potential moving decimations generating plaquette actions 
with large number of group characters and exhibit rapid approach to a unique renormalized 
trajectory. Fixing the critical coupling f3 c {N T ) at one value of temporal lattice length N T by 
MC simulation, the critical couplings for any other value of N T are then obtained by lattice 
matching of the block decimations. We obtain (3 C {N T ) values over the range N T = 3 — 32 and 
find agreement with MC simulation results to within a few percent in all cases. A similar 
procedure allows the calculation of string tensions with similarly good agreement with MC 
data. 
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1 Introduction 



In this paper we apply the renormalization group (RG) based technique of "lattice match- 
ing" via block transformations in pure lattice gauge theories. Instead of implementing block 
transformations by means of renormalization group Monte Carlo (MCRG) methods, how- 
ever, we employ explicit RG recursion relations of the "potential moving" type. The block 
transformations (decimations) implemented by these recursions are of course approximate 
but can, in principle, be systematically improved. They turn out to be surprisingly effective 
for various purposes. Here we use lattice matching of these decimations to obtain critical 
couplings and string tensions for the gauge groups SU(2) and SU(3). 

Specifically, we start with one critical coupling obtained by MC simulation at a certain 
temporal lattice size N T . We then employ lattice matching of our decimations to compute 
the critical couplings (3 C (N T ) at other values of iV T . We find that these values agree with 
the values obtained by MC to within at most a few percent. This procedure then affords a 
rather inexpensive method for obtaining critical couplings over a wide range of lattice sizes. 
In fact, since the recursions are locally specified, there appear to be no intrinsic lattice size 
limitations in the method. 

A different method for inexpensive estimation of critical couplings has recently been 
presented in [T]. In this approach the strong coupling expansion of the SU(N) theory is 
used to derive a 3-dimensional Z{N) effective Polyakov loop action. The easily ascertained 
critical values of the effective action parameters can then be mapped back to the critical 
couplings of the original SU(N) theory. It would be interesting to consider combining a RG 
decimation-based method such as the one presented here with that of p] to extend the range 
where the latter can be applied. 

Our lattice matching of decimations can also be applied, by a closely analogous procedure, 
to the computation of string tensions for various values of (3. Again, the computed values are 
in very good agreement with the MC data. 

The paper is arranged as follows. Our RG transformations and the resulting decimation 
recursions are formulated in section [2j The lattice matching method and its application via 
our decimations to the computation of critical couplings and string tensions are outlined in 
section [3l Our numerical results are presented in section [H Some concluding remarks are 
given in section 

2 RG decimations 

An RG block transformation with scale factor b maps a system on a lattice of spacing a to 
a system on a lattice of spacing ba. The flow in the space of interactions (couplings) under 
successive block transformations defines the RG flow of the system. 

Given the definition of some exact block transformation, its practical implementation, 
whether by analytical or numerical methods, generally involves some approximation or trun- 
cation. This is certainly the case if one aims at obtaining the effective action after each step. 
The standard method for implementation by numerical simulation is the MCRG method. 
Alternatively, one may incorporate some judicious approximations in the definition of the 
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transformation at the outset so that it becomes explicitly computable. Potential moving 
transformations are in this class. 

Partition a hypercubic (f-dimensional lattice of spacing a into <i-dimensional hypercubes 
of linear size ba (b an integer). The potential moving procedure [2] in the case of pure gauge 
theories consists of moving all plaquette interactions in the interior of each such hypercube to 
its (d — l)-dimensional boundary. After the move their strength is renormalized by some ad- 
justable factor to compensate for the move. Next the plaquettes inside this (d— l)-dimensional 
boundary are moved to its (d — 2)-dimensional boundary and similarly renormalized. Con- 
tinuing this process to its conclusion one ends up with the system on a lattice of spacing ba 
whose elementary plaquettes are tiled by at most b 2 plaquettes of the original lattice. 

This series of moves may be implemented in a number of somewhat different ways [3]. The 
simplest choice though is to first perform all moves in sequence and with common plaquette 
interaction renormalization £ after each move. This is the scheme we adopt here. It results 
in isotropic couplings in all directions, and plaquettes on the blocked lattice (of spacing ba) 
tiled by exactly b 2 plaquettes of the original lattice; and with each of these tiling plaquettes 
renormalized by a total factor £( d ~ 2 ). The integrations over those bond variables belonging 
to the tiling plaquettes, and interior to the blocked lattice elementary plaquettes, can now 
be performed and renormalized. This completes the block step yielding the theory on the 
blocked lattice. 

This procedure may be formulated as follows. Let, as usual, [/& £ G denote the bond 
variables, and U p = Ylbedp Ub their product around plaquette p. General elements of the 
gauge group G are denoted by U. Let A p (U p ,n) denote a plaquette action on lattice of 
spacing b n a and consider the character expansion 

exp (-A P (U P , n)) = ]T dj Fj(n) Xj (U p ) . (1) 
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The sum is over all inequivalent irreducible representations labeled by j, with Xj an d dj 
denoting the character and dimension, respectively, of the j-th. representation. From ([I]), 
using orthogonality of characters, one has 



Fj{n) = JdU exp (-A P (U, n)) 1 X *{U) , (2) 



where dU denotes Haar measure on G. In this paper we consider only G = SU (N),N = 2, 3. 
The action itself is of course completely specified by the set of Fj(n) coefficients and vice 
versa. It can be written in the general form 

A p (U p ,n) = ]T jty{n) ^IXjiUp) + XjiU" 1 )] (3) 
j j j 

with lj = 1 for self-conjugate and lj = 2 for non-self-conjugate representations. (For SU(2), 
in particular, L = 1 for all j.) In general we need consider actions of the form ([3]) with any 
(infinite) number of characters, i.e. couplings f3j. It is useful to define an effective coupling 
g( n ) characterizing a given action of the form ([3"]). With {t} denoting the SU(N) generators 
and \m\ = 1, we let 



oM ^iV _ d 2 A p {e^\n) 



ff (™) 2 d9 2 



(4) 
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(HJ is of course independent of the direction fa. In the perturbative regime this reduces 
to the usual definition of gauge coupling. In the non-perturbative regime any definition of a 
'coupling' is of course some scheme-dependent choice. We adopt (UJ) to track the RG evolution 
°f ©5 ©; it provides, in particular, a good parametrization of the renormalized trajectory 
(see below). 

A lattice block step b n a — > b n+1 a of the type described above can now be succinctly 
formulated as a prescription for obtaining the character expansion coefficients Fj{n + 1) in 
terms of the Fj{n): 



FAn + 1) 



dU 



J2 dtFiWxiiU) 



d-2) 



(5) 



The inner bracket factor results from the symmetric potential moves described above. Subse- 
quent boundary integrations modify each resulting expansion coefficient by a further amount 
controlled by the parameter r. Specification of £ and r completes the block step. 

Our decimation transformation b n a — > b n+l a is defined by ([5|) with [3]: 



c 



1 - eg 



(n)2 



b 1 - eg 



(n)2 



(6) 

(7) 



where c is an adjustable decimation parameter. For sufficiently large effective coupling values 
the </ n ) 2 dependence in © - (J7J) has to be suitably leveled off [7], since we require £ > 0, 
r > 0. This is, however, not explicitly indicated here as this regime is well outside the range 
of couplings encountered in our applications of the recursions ([5]) below. 

In the following c, for given b, will be treated as a decimation parameter tuned for op- 
timization of the procedure outlined in the next section. It is interesting to note, however, 
that this turns out to give values (section HJ) not far from the weak coupling computed values 

[HIS]. 

It is convenient to work with normalized coefficients fj = Fj/Fq by factoring out the 
trivial representation coefficient in ([1]); the correspondingly normalized action differs by a 
trivial shift of the constant (trivial character) part in ([3]). Effective couplings (j4|) are also 
conveniently computed directly in terms of the {fj}. 



3 Lattice matching of decimations 

Given a ci-dimensional lattice system with action A(K) defined by a set of couplings K = {Ki} 
RG block transformations by a scale factor b generate a flow in action space: 

K (o) K W ^R 2 ^ ► (8) 

where = {K^} denotes the couplings after n blocking steps from the initial point 

Since the physical correlation length remains of course fixed, the (dimensionless) 
lattice correlation length £( n ) at step n gets rescaled as £( n ) = jb n . 
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The resulting RG flow is toward a fixed point along irrelevant directions (couplings) and 
away from the fixed point along relevant directions (couplings). Irrespective of the starting 
point then, the flow, after a sufficient number of blocking steps, will approach the unique 
Wilsonian renormalized trajectory (RT) emanating from the fixed point along the relevant 
directions. 

Consider two sets of couplings K and K' . If the two RG flows starting from them end up 
at the same point on the RT after the same number of blocking steps n, then, since at the end 
point £( n ) = £'( n \ the correlation lengths = £ and = £' at K and K', respectively, 
must be equal; and since the physical correlation length is constant, K and K' must also 
have the same lattice spacing a. By the same token, if the flows from K and K' reach the 
same point on the RT after n and (n — m) steps, respectively, then the lattice correlation 
lengths at K and K' must be related as 

£ = b- m t- (9) 

and the lattice spacings a and a' as 

a' = b m a. (10) 

To identify such pairs of couplings we need to ascertain that, after n and (n — m) RG 
steps respectively, the same point is reached on the RT. This can be done in two ways. One 
is to show that A(K^) = A(K'( n ~ m ^). This requires that one obtain the blocked action at 
each step. Another way is to show that the expectations of every operator, measured after 
performing the corresponding number of blocking steps from the initial two actions, agree. 
Either way, blocking n times from a starting point K, and then adjusting another starting 
point K' so that after (n — 1), or, more generally, (n — m) times, matching is achieved is 
referred to as two-lattice matching [8]. 

If blockings are performed numerically by MCRG, the second method appears easier to 
use. Obtaining the blocked action can be difficult, whereas it is possible, at least in principle, 
to generate a Boltzmann-weighted configuration ensemble for the blocked action by instead 
blocking the configurations of an ensemble generated from the original action. These can 
then be used to measure observables [9]. In practice, of course, due to lattice size limitations, 
only a rather small number of block steps is possible by MCRG, so getting close enough to 
the RT is not guaranteed. In this connection, since the location of the fixed point is block 
definition dependent, appropriate fine-tuning of free parameters in the block transformation 
definition can be crucial for achieving rapid approach in few steps. 

In this paper we employ two-lattice matching with RG block transformations implemented 
by the recursions © described above. They can be explicitly evaluated to any desired 
accuracy on lattices of any size, so no inherent limitations due to finite size arise. The blocked 
action resulting after each RG step is explicitly obtained, so it can used to ascertain approach 
to the RT and perform two-lattice matching. The transformations contain one parameter (cf. 
& - ©); which, as already mentioned, is fixed for optimization of the matching. 

In the following the starting action (n = 0) will always be taken to be the fundamental 
representation Wilson action. Other choices such as mixed actions containing several repre- 
sentations can be treated in exactly the same way. The flow under successive decimations 
reaches a unique RT irrespective of such a choice, though of course the number of steps needed 
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to reach it depends on the initial point in action space. The important feature characterizing 
these decimations is that, regardless of the choice of the initial plaquette action, a single step 
suffices to generate an action of the form (|3j) generally containing the full (infinite) set of 
representations. Flow in such an infinite-dimensional interaction space makes it possible to 
avoid getting stuck at (finite-dimensional) lattice artifact boundaries. 

With the fundamental representation Wilson action as the starting point we find that the 
approach to the unique RT is very rapid; it generally takes only two steps to get to it. This 
is illustrated in Figured) 




Figure 1: RG flow and rapid approach to a unique renormalized trajectory starting from the 
SU(2) fundamental representation Wilson action with f3 = 4 (green dots) and (3 = 2.5 (red 
dots). The first three non-trivial (normalized) expansion coefficients in ([1]) are shown. 

A good way to parametrize points along the RT is by the effective coupling (J5|) of the 
action corresponding to each RT point. If, starting from some Wilson action coupling /3, 
after n steps the point /3^ n '(/3) lies on the RT, subsequent RG steps generate a sequence of 
points f3^ n+l \ f3^ n+2 \ ■ ■ ■ hopping along the RT. With scale factor 6 = 2, and for all large and 
intermediate values of /3^ n \ the effective beta function is varying slowly enough for a linear 
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interpolation to provide an excellent approximation to the RT points lying between pairs of 
neighboring points P^ n \(3^ n+1 \ So we write 

p(n+h) = p{n) + h (p(n+l) _ p(n)^ ^ < 7i < 1 . (11) 

This may be viewed as performing n steps followed by a step with fractional scale factor to 
reach a point j3^ n+h ^ lying between point and /3( n+1 ) on the RT. (The block transformation 
rules ([5j)-((TJ) can indeed be formally extended to non-integer scale factor, but we need not 
make actual use of this here.) (jlip asserts that the location of this intermediate point is 
accurately given by linear interpolation. 

Though explicitly computable to any accuracy, our decimations are of course approximate. 
They become exact in the strict /3 — > oo limit. Computation of the step scaling function 
(beta function) from the decimations in the weak coupling scaling region reproduces the 
perturbation theory prediction to within 2% — 3%. The next question to be probed by 
further computation then is how close an approximation these decimations give to the exact 
RT in the transition to the non-perturbative regime. MCRG construction of blocked actions 
[10] shows that one-plaquette terms with a large number of characters are the most relevant 
action terms for long-scale dynamics. This is precisely the type of action resulting from our 
decimations. 

In the following two-lattice matching of our RG decimations is used to obtain critical 
couplings and string tensions for the SU(2) and SU(3) gauge theories. 



3.1 Critical couplings 

Consider the (3 + l)-dimensional lattice theory at physical temperature T. Since T = l/aN T 
for lattice of time extent N T and spacing a, lattice of extent N^, spacing a' is related by 

a ' = w a - (12) 

T 

If after blocking the two lattices n and n' times, respectively, the two flows reach the same 
point on the RT trajectory, (fT0|) and (fl~2j) imply 

'=•*(§)■ o*> 



n — n 



If, in particular, T = T c , one has 

^(/3 c (N T ))=^(/3 c (N! r )). (14) 

( I13p and (|14|) afford a straightforward way of obtaining critical couplings by matching 
once one such coupling is known. Assuming f3 c {N T ) known, it is convenient to simply choose 

n = log ft N T 

n' = lo gf X (15) 

so that (fT3l) is satisfied, n, n' must be large enough to be on the RT, but this is not a problem 
since one can always simply add a common integer to the r.h.s. of both equations in ()15p . 
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Also, note that the so-chosen n or/and n' may turn out to be non- integer. In such a case, one 
performs [n] and [n] + 1 steps, where [n] is the nearest integer to the chosen n from below, 
and uses (jlip - and similarly for n' . With n, n' and /3 C (-/V T ) given, (|14p is then solved for 
j3 c (N' T ), i.e., the starting point of the flow on the N' T lattice is adjusted to satisfy (fT4|) . 



3.2 String tensions 

The string tension <r(/3) = a 2 a is another quantity that can be obtained at different couplings 
by the same method. Consider two RG flows with starting point the Wilson action at 
couplings /3o and j3\ and ending up at the same point on the RT after no and n\ steps, 
respectively. Then 

/? (no) (/?o) = /3 (ni) (/3i) (16) 

and 

ai ^ = b^°- n ^aoV^ (17) 

by CED. 

Suppose we know <5"(/?o) • Choose no large enough to be on the RT. Then n\ is fixed so 
that (|16|) is satisfied. In practice this is done by performing n and n + 1 decimation RG steps 
with initial coupling Pi, where n is determined so that 

(p^iPo) - /3 (n+1) (/3i)) (/5 (no) (/3o) - /3 (n) (/3i)) < . (18) 
Then, by (lllj) . n\ satisfying ()16f) is given by n\ = n + h with 

i /3 ( " o) (/3o) - /? (ri) (/3i) 



/3(«+ 1 )(/3 1 )-/3W(/3 1 ) 
y/&(/3i) = a\yfa is then obtained directly from ([1] 



(19) 



4 Results for critical couplings and string tensions 

For SU(2) we typically use fifty group characters in the expansions (pQ). This implies for, say, 
(3 = 5 omitted higher character coefficients fj = Fj/Fo, and accompanying bounds on the 
series remainder, of the order of 10 -45 . For SU(3) we truncate ([T]) at characters j = (p,q) 
with p > 20, q > 20; this implies remainders at /3 = 10 of less than 10~ 12 . Iteration under 
([!]) results into decreasing expansion coefficients. Errors due to truncation in the character 
expansions (pQ) are thus totally negligible. 

The scale factor is always taken to be b = 2. The only adjustable parameter in the 
decimation recursions ([5]) - ([7J is c which is tuned for optimized matching. We set c = 0.10 
in the case of SU(2) and c = 0.24 in the case of SU(3). With no other parameters present, 
straightforward numerical evaluation of the recursion relations can then be implemented. 

We take one value of (3 C (N T ) from MC data, which serves to fix the scale. We then apply 
the procedure of section I3TT1 to obtain critical coupling values for other lattices. Results for 
SU(2) are shown in Table [TJ Two sets of computed /3 C values are shown in Table [1] (columns 



8 



N 


R 
Pc 


R 
He 


R (MC) 


Q 

o 




z.lyo i 


z.i f Oo^oU J 


A 
*± 


9 9QDQ 


9 9QQ1 


9 9QQ1 (C\9\ 


5 


2 3600 


2 3683 


2 3726(45'! 


6 


2.4175 


2.4258 


2.4265(30) 


8 


2.5097 


2.5180 


2.5104(02) 


12 


2.6355 


2.6440 


2.6355(10) 


16 


2.7275 


2.7361 


2.7310(20) 


32 


2.9487 


2.9574 





Table 1: Critical couplings (3 C (N T ) for SU(2) computed from lattice matching of decimations. 
Column 1 and 2 show the values obtained for two different choices (underlined entries) of 
the one data point taken from MC data (see text). Column 3 shows the values from MC 
simulations for comparison. 





/3 C 


/5 C 


ft(MC) 


4 


5.6501 


5.6329 


5.6925(002) 


6 


5.8941 


5.8773 


5.8941(005) 


8 


6.0773 


6.0595 


6.0010(250),6.0625(18) 


10 


6.2018 


6.1837 


6.1600(70) 


12 


6.3084 


6.2900 


6.2680(120),6.3385(55) 


14 


6.4015 


6.3830 


6.3830(100) 


16 


6.4845 


6.4658 


6.4500(500) 


32 


6.9024 


6.8829 





Table 2: Critical couplings f3 c (N T ) for SU(3) computed from lattice matching of decimations 
and comparison with MC simulation data. Same format as in Table HJ 

1 and 2) corresponding to two different choices of the MC data point (underlined entries). 
The table also shows comparison with the values obtained by MC simulation [11] - [13] , [1] in 
each case (column 3). The agreement is remarkably good - typically of the order of 1% — 2%. 

Results for critical couplings in the SU (3) gauge theory are displayed in Table [2j Agree- 
ment with MC simulation data [11] , [1] is again very good, typically within a few percent. 

String tensions in SU(2) obtained by the method of section [3T21 are displayed in the same 
format in Table Again, two sets of values are shown (columns 1 and 2) corresponding to 
two different choices (underlined entries) of the MC data point used as initial input. The 
corresponding results in the case of SU (3) are shown in Table HI Good agreement with MC 
data [13] - [15] is again obtained in all cases. 

5 Conclusions 

The RG decimation recursion relations given in section [3] were used in conjunction with two- 
lattice matching to compute critical couplings and string tensions in SU(2) and SU(3) pure 
lattice gauge theories. The decimations contain only one adjustable parameter that was fixed, 
in the case of each group, to an optimized value given in section [H Using one initial value 
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ft 
p 






n frr(MC\ 


9 9 


u.ouiy 


U.01D1 


u.^oyu^iuu ) 


9 ^ 

z.o 


U.OO04 


U.o / OD 


u.ooyu^ou ) 


9 /I 


n 9R1 Q 

u.zoiy 


u.zoyo 


U.ZOOU^ZU ) 


2 5 


f) 1903 


1 957 


1Q05(08"l 


2.5115 


0.1836 


0.1888 


0.1836(13) 


2.6 


0.1373 


0.1415 


0.1360(40) 


2.7 


0.1002 


0.1031 


0.1015(10) 


2.74 


0.0884 


0.0911 


0.0911(08) 


2.85 


0.0622 


0.0641 


0.0630(30) 



Table 3: String tensions a^fo for SU(2) computed from lattice matching of decimations. 
Column 1 and 2 show the values obtained for two different choices (underlined entries) of 
the one data point taken from MC data (see text). Column 3 shows the values from MC 
simulations for comparison. 



(3 


a^fa 


ay/a 


ay/a(MC) 


5.54 


0.5580 


0.5878 


0.5727(52) 


5.6 


0.5070 


0.5295 


0.5295(09), 0.5064(28) 


5.7 


0.4205 


0.4264 


0.4099(12), 0.3879(39) 


5.8 


0.3486 


0.3508 


0.3302(15) 


5.9 


0.2919 


0.2931 


0.2702(19) 


6.0 


0.2465 


0.2433 


0.2269(62), 0.2209(23) 


6.2 


0.1698 


0.1671 


0.1619(19), 0.1604(11) 


6.4 


0.1214 


0.1180 


0.1214(12), 0.1218(28) 


6.5 


0.1010 


0.0983 


0.1068(09) 


6.8 


0.0616 


0.0599 


0.0738(20) 



Table 4: String tensions a\fa for SU(3) computed from lattice matching of decimations. 
Same format as in Table El 

obtained by MC simulation, critical couplings and string tensions were then obtained for a 
variety of other lattices by lattice matching of our decimations. The results were found to 
be in very good agreement with those obtained by MC simulation. The method evidently 
provides a cheap way of quickly obtaining accurate predictions for these quantities for a wide 
range of lattice sizes. 

Critical couplings and string tensions are quantities pertaining to long-distance non- 
perturbative dynamics. The actions evolving under the decimations are plaquette actions 
with a large (infinite) number of representations. As mentioned above, MCRG constructions 
of blocked actions [10J indicate that these are the action terms most relevant for long distance 
dynamics. This may be one reason underlying the method's apparent efficacy. 

There are two directions in which this work could be further pursued. One is to consider 
more general block transformations. The decimations employed here may indeed be viewed 
as special cases of more elaborate blocking schemes. These will, in general, involve additional 
decimation parameters, but are likely necessary for computation of observables over differ- 
ent length scales. The other direction is the inclusion of fermions. Block transformations 
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involving fermions present a generally much harder problem. Use of relatively simple block 
schemes in a somewhat different context, however, have given some early encouraging results 
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